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We study the phase diagram of the frustrated Heisenberg model on the triangular lattice with 
nearest and next-nearest neighbor spin exchange coupling, on 3-leg ladders. Using the density- 
matrix renormalization-group method, we obtain the complete phase diagram of the model, which 
includes quasi-long-range 120° and columnar order, and a Majumdar-Ghosh phase with short-ranged 
correlations. All these phases are non-chiral and planar. We also identify the nature of phase 
transitions. 


I. INTRODUCTION 

Quantum magnetism in reduced dimensions gives rise 
to a fascinating range of behaviors.^^^In one-dimensional 
(ID) systems, powerful analytical and numerical methods 
have allowed a deep understanding of phenomena such 
as fractionalization,!^ di mer ization,^ and symmetry pro¬ 
tected topological order.l^In two dimensions (2D) there 
remain many m ore open problems such as un derst and- 
ing spin liquids,intrinsic topological order, ESEI] and 
the connection between exotic mag netic phases and un¬ 
conventional superconductivity.*^*^ Few-leg ladders are a 
vital intermediate class as they allow for the application 
of accurate numerical methodJ^^ available for large ID 
systems, while also providing important insights int o new 
physics occurring in the crossover to two dimensions !^^ * ^^ * 

In a unfrustrated system, such as the nearest neigh¬ 
bor Heisenberg model on the square lattice, all of the 
terms in the Hamiltonian can be minimized simultane¬ 
ously. This tends to favor long-range order in 2D. There¬ 
fore, frustrated systems are excellent candidates in which 
to searc h for exotic phases of matter without conventional 
ordering.*^ 

The spin-1/2 triangular Heisenberg model (THM) 
is a prototypical model for frustrated magnets in 
two dimension.!^ In 1973, Andersorfi^ suggested that 
the resonating-valence-bond (RVB) state could play 
a pivotal role in the description of novel magnetic 
materials, and his conjecture that ground-state of the 
spin-1/2 THM would be an RVB state provoked much 
interest. However studies of this model hav e failed 
to find an RVB state and the evidenc^^^**^ is now 
very strong that for the pure isotropic model with 
nearest-neighbor (NN) interactions, the ground-state is 
a 120° magnetically ordered state.*^* Variants of the THM 
describe some properties of organic materiald^ such 
as k-(BEDT-TTF) 2 Cu 2 (CN) 3 , EtM e.3Sb[ Pd(dmit)2]2, 
EtMe3P[Pd(dmit)2]2, Mo3S7(dmit) J^° l ^d an d al so 
some quasi-two-dimensional inorganic material j^ * ^ * ^^* *^ 
such as RbFe(Mo 04 ) 2 , Ba 3 CoSb 209 , Cs 2 CuBr 4 , and 
CS 2 CUCI 4 . 

In ID the prototypical frustrated system is the zig-zag 
chain, which has an exact solution at the Majumdar- 


Ghosh point ,1^ characterized by long-range dimer order 
and a two-fold degenerate ground-state. As we show be¬ 
low, an NNN Majumdar-Ghosh phase is stabilized in a 
large region of the phase diagram of the 3-leg triangular 
ladder. 

So far the THM has been mostly considered with only 
NN exchange coupling, but additional interactions or 
anisotropies may stabilize exotic states. A natural choice 
for an additional interaction, while retaining isotropy, 
is a next-nearest neighbor (NNN) coupling to add fur¬ 
ther frustration effects. In this paper we study the J1-J2 
THM on a width 3 cylinder as a simplified version of the 
full 2D model, but readily accessible to numerical meth¬ 
ods. The ladder model has clear connections to the 2D 
THM and also extrapolates smoothly to the Majumdar- 
Ghosh point of the zig-zag chain. The J1-J2 THM in 
2D has been previously studied using semi-clas sical spi n- 
wave theories (SWT) and exact diagonalization,*ISlM22l29| 
but these studies did not cover the physics of the whole 
phase diagram. Recently, a coupled cl uster studji^ and 
quantum Monte Garlo (QMG) studie^^i*^ have identi¬ 
fied a phase in this model that is a candidate for a spin 
liquid. Magnetically ordered states in a variety of clas¬ 
sical 0(3) models with J 1 -J 2 -J 3 interactions have been 
studied in [331 finding several different ‘regular magnetic 
orders’ relevant to the triangular lattice, including planar 
and non-planar 120° states. 

The lattice we consider is shown in Pig. The Hamil¬ 
tonian is 

h = j,J2s,.s, + j, S-S.-. W 

(hi) ((hi)) 

where {i,j) (((*, j))) indicates that the sum is over all NN 
(NNN) couplings. To cover the full range of couplings, 
we introduce the following parameterisation: 

Ji = J cos 9, J2 = J sin 9 . (2) 

Where the J is the unit of energy and henceforth we 
fix J = 1.0. The main difficulty in studying the lattice 
shown in Pig. [^is frustration. The lowest energy state of 
the AFM Heisenberg model on a square lattice has Neel 
order.*^ This cannot be formed on an equilateral triangu¬ 
lar lattice, and as a result there is competition between 
















2 



FIG. 1. (Color online) Lattice structure and mapping of the 
3-leg cylinder to the one-dimensional chain employed by the 
finite DMRG calculations. Spins sit on vertices. The lattice 
is tripartite and the sub-lattices are labeled A, B, and C. ai 
and a 2 are lattice vectors. 

terms in the Hamiltonian, Eq. ([^, and they cannot si¬ 
multaneously minimize their local energy. Therefore it is 
clear that the lowest energy state must be a compromise, 
such as the 120° state. The 120° state on the triangular 
lattice is less stable than the Neel state on the square 
lattice,^ as the sublattice magnetization of the triangu¬ 
lar lattice is significantly reduced compared to its classi¬ 
cal value. Because of this reduced stability inherent to 
the triangular lattice, upon perturbing the Hamiltonian 
one may expect to see a variety of new phases. 

There have been several numerical studi es of the T HM 
in the past. Exact diagonalization methodJi^IllHllIH 
fer from exponentially growing size of the Hilbert space, 
which is especially a problem in two or more dimen¬ 
sions, while QMC technique^^ suffer from the sign prob¬ 
lem for frustr ated l attices, and projected entangled pair 
states (PEPSjpS*^ for this model is complex and com¬ 
putationally costly, although in principle PEPS has good 
computational scaling properties in 2D. More recently 
some numerical methods have been developed that are 
especially useful for frustrated systems and applied to 
the THM. For example, large-scale parallel tempering 
Monte CarlcP^ and some tensor networks methods in¬ 
cluding entangled-plaquette stated and multi-scale en¬ 
tanglement renormalization ansatz (MERA).I^ 

On the other hand, matrix product states (MPSjP^ 
have been around in various guises for a long time and are 
a good representation of the ground-state of ID and few- 
leg ladders. MPS exploits the locality of the interactions 
for ID ladders, and computes a truncated Hilbert space 
that is well-suited for describing ground-states, as it sat¬ 
isfies the area law for the bipartite entanglement (see [36] 
and references therein). In particular, the density matrix 
renormalization groupS"^ (DMRG) method for Hnding 
the variational ground-state is mature and highly effi¬ 
cient. A recent studjS^ of different numerical methods 
suggested that two-dimensional DMRG could be “one of 
the most powerful methods” for studying quantum lat¬ 
tice systems. 

A study using same method as this paper derived 


the phase diagram of the J 1 -J 2 Heisenberg model on 
a Kagome lattice, which also contains a rich variety of 
phases,!^ including a spin liquid and magnetically or¬ 
dered states. The THM on a 3-leg ladder has been pre¬ 
viously studied for anisotropic NN interactions (J 2 = 0) 
with a magnetic held,S^ and a phase diagram has been 
obtained. At the isotropic point, corresponding to d = 0 
in our notation, it was shown that the introduction of 
a magnetic held, causes the 120° state to 

evolve into “commensurate planar” phases with Y- and 
V-shape spin ordering on either side of a 1/3 magnetiza¬ 
tion plateau.l^Sl 


II. METHODS: MPS AND DMRG 

In this paper we employ the MPS ansatz, keeping up to 
m = 1000 basis states, using the DMRG method for ob¬ 
taining the ground-state wavefunction. The Hamiltonian 
has SU{2) symmetry, 

[H,S] = 0. (3) 

Exploiting this symmetry in the calculations gives a 
signihcant improvement in efficiency, by reducing the 
dimension of the computational Hilbert space. Using 
m = 1000 SU(2)-symmetric basis states is equivalent 
to m « 3000 states with no (or just Abelian 17(1)) 
symmetry. We performed both hnite DMRG and inh- 
nite DMRG)^ (iDMRG) calculations. The latter exploits 
translational symmetry available in the thermodynamic 
limit. 

Because MPS is fundamentally a ID ansatz, to apply 
it to 2D models a mapping is necessary. We map the 
ID chain of spins into a size N = L x 3 chain as shown 
in Fig. where L is the length of the 3-leg cylinder. 
The computational cost will scale approximately linearly 
with length, but still exponential with width, which is a 
limitation of this method. The model is on a cylinder, 
i.e. we use open boundary conditions (OBG) in the long 
(horizontal) direction and periodic boundary condition 
(PBC) in the short (vertical) direction. Using this wrap¬ 
ping, the lattice is conventionally a YCLx3 and there 
are {3N — 12)-bonds available for each of NN and NNN 
couplings. From finite size scalings of the energy and 
other order parameters, we found that sizes up to 30 x 3 
are large enough to scale finite results properly into the 
thermodynamic limit. 

In the case of iDMRG, one can classify all possible 
wrappings of the triangular lattice on an infinite 3-leg 
cylinder, using a standard notation developed for single¬ 
wall carbon nanotubes.^ We use the wrapping vector 
Cg = (—3, 3) in this notation. The unit vectors, ai and 
a 2 , used to specify Cg, are shown in Fig.j^ Cg preserves 
the tripartite symmetry on the infinite lattice. The pitch 
angle of this wrapping method is t/g = 90°. The matrix 
product operator^ (MPO) representation of the Hamil¬ 
tonian has a 3-site unit-cell in the direction of Cg. One 
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FIG. 2. (Color online) Energy versus variance for a 30 x 3 
cylinder with 6 = 25°. The linear extrapolation gives a good 
approximation for the exact ground-state energy. 



can show that Cq is the shortest possible wrapping vec¬ 
tor that preserves tripartite symmetry. The minimum 
unit-cell of the wavefunction however is 18 sites, as the 
smallest even size that preserves tripartite symmetry. 


FIG. 3. (Color online) The calculated phase diagram of the 
J 1 -J 2 THM on a 3-leg cylinder. The phase transitions are 
indicated to a resolution of 0.5°. All transitions are second- 
order except for 0 = — which is first-order (marked by a 
thick black line). 


A. Error Analysis 

We use the variance to calculate systematic errors in 
the DMRG results. E.g. in the case of energy, we have 
(t|; = {'4>v\{H — EY\ijjy). For energy errors, one needs 
to plot energy versus variance step-by-step for different 
numbers of states, m. The behavior of E versus is 
expected to be linear. Any significant deviation from 
this linearity indicates that the DMRG calculation has 
not converged, possibly due to an insufficient number of 
basis states. An example of this calculation is shown in 
Fig.d for a 30 X 3 lattice with 6 = 25°, for m between 
500 and 1000. The ground-state energy extrapolated to 
the m —)■ 00 limit is i?o[oo] = —46.94877331266(2). This 
method is similar to, but more robust than, the energy 
versus truncation error scaling that is typically used in 
DMRG calculations.!^ 

Throughout this paper the results are all converged 
with relative errors ^ 10“^^ — 10“®. Errors are smaller 
than symbol size for all plots except for the finite-size 
extrapolation of the spin gap in Fig. |15[ a), where we 
show the error bars. 


III. PHASE DIAGRAM 

The calculated phase diagram of the J 1 -J 2 THM is 
shown in Fig. The dominant short-range ordering is 
sketched in the form of triangular or rhombic plaquettes. 
The model contains four well-defined phases. The dif¬ 


ferent phases were determined by studying the ground- 
state energy (Sec. IIIB), spin-spin correlation functions 


(Sec. |IV[ ), the chirality (Sec. 120° order parameter 
(Sec. |VI| , spin gap, dimer order parameter, and Binder 
cumulant (Sec. VII). 

In order to better visualize the nature of the short- 
range correlations in each phase. Fig. shows the NN 
and NNN bonds, colored according to the value and sign 
of the spin-spin correlation. The four phases are: 

1) A 120° state (cf. Fig. [4|(a)) that exists in 
the fourth quadrant of Fig. The 120° state 
is critical (see Sec. VII below), with power-law 
correlations and gapless excitations. This is 
in contrast to the 2D model, which has long- 
range tripartite magnetic ordering. However long- 
range magnetic ordering is forbidden in our 3-leg 
cylinder due to the Mermin-Wagner theorenPS, 
which excludes SU (2) symmetry-broken long-range 
order in ID. The 120° state is Ge rotational 
symmetric, parity-symmetric (P-symmetric), time- 
reversal-symmetric (T-symmetric), and planar (see 
Sec.j^below). This state persists in the first quad¬ 
rant up to a quantum critical point at non-trivial 
9c — 6.5°. The existence of 120° state is consis¬ 
tent with spin-wave results of Jolicouer et al.!^^ Al¬ 
though the transition point of spin-wave calcula¬ 
tions is located at 9swt = tan“^ | ~ 7.125° com¬ 
pared to our value of 6.5°. 

2) Upon increasing 9 > 6.5°, the system changes 
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FIG. 4. (Color online) Visualizations representative of differ¬ 
ent phases of the THM on a cylinder. Thicker lines represent 
stronger bonds with red indicating antiferromagnetic bonds, 
blue indicating ferromagnetic. 

phase to a two-sublattice commensurate spin state 
with a columnar structure (cf. Fig. |^b)), which 
is also gapless. This is in contrast to the 2D 
model, which has long-range columnar order.l^SHUl 
However this ordering is forbidden in ID for the 
same reason as the 120° state. The columnar state 
is quasi-long-range, Ce rotational symmetry bro¬ 
ken, P-symmetric, T-symmetric, and planar. This 
phase can be thought of as a planar version of the 
standard G-type antiferromagnetism.l^ 

3) At 0 ~ 70.0°, there is a phase transition to a 
NNN Majumdar-Ghosh state. In this phase 
the system forms strong AFM bonds (dimers) from 
the J 2 interaction, in a way that only one of the 
two available double-counted edges of any sub¬ 
lattice triangular plaquette has a strong (singlet) 
bond, leading to a two-fold degenerate ground-state 
with long-range dimer order. The phase becomes 
a perfect Majumdar-Ghosh singlet-dimer state at 
exactly 0 = | (cf. Fig. |^c)). We will see be¬ 
low that at this point the model is composed of 
three uncoupled and dimerized sublattices in the 
form of 2-leg spin ladders, which is a direct con¬ 
sequence of the PEG in the short direction of the 
lattice. The Majumdar-Ghosh state is robust to 
small perturbations when one turns on the Ji in¬ 
teractions, and evolves into the general form shown 
in Fig. |4];d), with weak NN bonds, either antifer¬ 
romagnetic or ferromagnetic corresponding to the 
sign of Ji. The Majumdar-Ghosh state has short- 
ranged correlations (cf. Fig.|^, and is Ge rotational 
symmetry broken, P-symmetric, T-symmetric, and 
planar. The state persists throughout a large re¬ 


gion in the second quadrant of the Fig. 


4) Upon further increasing of 6 , the system undergoes 
a second-order phase transition at 9c = 152.0° (see 
Sec. VII below). In a narrow region, 152° < 9 < 
165° of Fig. 1^ the ground-state is a partially polar¬ 
ized ferromagnet that saturates to complete fer¬ 
romagnetism for 9 > 165°. 


A. Limiting Cases 


In our parameterization of the Hamiltonian, 0 = 0 is 
equivalent to J 2 = 0, and is simply the nearest-neighbor 
model. The ground-state is the 120° state, in agreement 
with the semi-classical approach,li^with wave vector Q = 
{2tt/V3, 27r/3) in our notation. 

For 9 = 90° (Ji = 0), the model has only NNN in¬ 
teractions. This state is composed of three uncoupled 
spin ladders, one in each tripartite sublattice, forming 
a perfect Majumdar-Ghosh state of alternating singlet 
dimers.l^ The formation of this phase is a direct con¬ 
sequence of the 3-leg form of the lattice. Fig. which 
is wrapped around a cylinder resulting in 3 independent 
zig-zag spin chains with NN coupling J 2 , and the double- 
counted bonds around the periodic boundary give a NNN 
coupling of 2 J 2 . Thus, this state appears because of the 
restricted geometry of the 3-leg ladder. On the other 
hand, in the 2D limit the Hamiltonian is instead three 
copies of the 9 = 0 model, hence the ground-state will 
contain three copies of the 120° state, one on each sub¬ 
lattice, and a small Ji will couple the otherwise indepen¬ 
dent sublattices. Thus the small-Ji behavior for few-leg 
ladders is rather different to the bulk 2D behavior. 

For 0 < 0 < tan“^(|) the 2D model at the classi¬ 
cal level {S —>■ 00 ) has a 120° ground-stat^^ and for 
9 > tan“^(|) it has a 4-sublattice AFM Neel phase with 
an infinite manifold of degenerate ground-states, selected 
by the “order from disorder” mechanism. Quantum fluc¬ 
tuations break this degeneracy, and the quantum model 
has a two-sublattice columnar (collinear) Neel state.l^^^ 
It is worth mentioning that the selection of the collinear 
order from the 4-sublattice classical order, could be un¬ 
derstood analytically using group-symmetry analysis 

For the classical J 1 -J 2 THM, if one enforces the tripar¬ 
tite symmetry everywhere using a repeated 3-site unit¬ 
cell, the ground-state phase is simply ferromagnetic (FM) 
for Ji < 0 and the 120° state for Ji > 0, independent of 

J2- 


B. Ground-State Energy 

In this section, we benchmark our results for the 
ground-state energy per nearest-neighbor (Ji) bond, Eq- 
This is shown in Fig. The energy per NN bond in the 
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FIG. 5. (Color online) (a) Ground-state energy of the THM 
on a 30 X 3 cylinder. The transition at 0 = — ^ is sharp and 
consistent with a hrst-order phase transition (cf. Sec. VIII. 
The transition near 6 « 160° snggests a second order phase 
transition (cf. Sec. |VIl| ). Brown stripes indicate phase transi¬ 
tions. (b) Comparison of the energies of Lanczos and DMRG 
for the THM on a 4x3 lattice. Lanczos results are from Joli- 
coeur et al.^DMRG and Lanczos show excellent agreement, 
but the boundary conditions clearly have a signihcant effect 
on the energy of this small lattice. 



e 


FIG. 6. (Color online) Average of squared magnetization 
per plaquette, Eq. on a 30x3 cylinder. The turquoise 
stripes denotes the partially-polarized region at the second- 
order transition into the ferromagnetic state. 


FM is second-order, which we verified by calculating the 
magnetization (see Sec. VII below). 

Fig.j^b) is a comparison of DMRG energies with Lanc¬ 
zos results of Jolicoeur et al.^ They simulated the same 
model on a 12-site lattice with PBC in both directions, 
which is equivalent to a 4 x 3 torus in our representation. 
The choice of wrapping vector around the torus has lit¬ 
tle effect as long as lattice translational and tripartite 
symmetry are preserved. Our DMRG results are in very 
good agreement with these Lanczos results. 

Table |l] is a comparison between our DMRG energy 
and results from previous calculations for 0 = 0, i.e. the 
NN model in the 120° phase. For this point we performed 
a larger size calculation on a 60 x 3 cylinder, as there is 
no NNN frustration and the DMRG is easier to converge. 
The results in Table |T] suggest that the THM on a cylinder 
is a good approximation for the full 2D model. 


C. Local magnetization 


fully polarised ferromagnet is, 

^ (sin 0-I-cos 0) , (4) 

which is shown in turquoise in Fig.j^a). 

There is a sharp transition appearing at 0 = —7r/2, 
coinciding with the change from FM to 120° state. The 
cusp suggests a first-order phase transition, which is con¬ 
firmed by the local magnetization and order parameters. 
This is the only first-order transition that we find in the 
model, and is indicated by the thick black line in Fig. 
On the right-hand side of Fig. [§a), in the vicinity of 
0c = 152.0°, the derivative is continuous indicating that 
the transition from the Majumdar-Ghosh state to the 


The squared magnetization per plaquette, is pre¬ 
sented in Fig.|^ This is calculated from the square of the 
local magnetization on a single plaquette, 

= ^ E (SA + SB + Sc)^ (5) 

where the sum is over all Np plaquettes with vertices 
A, B, and C from their respective sublattice. The 
turquoise lines in Fig. indicate the region where we 
found a partially-polarized ferromagnetic ground-state. 
The rapid but smooth change in local magnetization in 
this region is consistent with a second-order phase tran¬ 
sition. 
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TABLE 1. Comparison of the ground-state energies from different methods for the THM at 0 = 0. 


Method 

Gluster size 

Boundary conditions 

Energy per bond, e® ° 

DMRG (this work) 
DMRG (this work) 

DMRG (this work) 

4x3 

60x3 

extrapolated to thermodynamic limit 
(L X 3 lattices with L < 60) 

torus 

cylinder 

cylinders 

< -0.20164623520324(1) 
< -0.19053054(3) 

-0.189(2) 

iDMRG (this work) 

inhnite 

3-leg cylinder 

< -0.189715084187(2) 


Ising expansiorP^ 

Entangled plaquette state^^ 

Coupled cluster methocP^ 

Numerical diagonalizatiorP^ 

QMCP 

SWlE^ 

MERAESI 


iV = 12 

theoretically thermodynamic limit 
extrapolated to thermodynamic limit 
(clusters up to size N = 324) 
extrapolated to thermodynamic limit 
extrapolated to thermodynamic limit 
(clusters up to size N = 27) 
extrapolated to thermodynamic limit 
theoretically thermodynamic limit 
extrapolated to thermodynamic limit 
(clusters up to size N — 114) 


torus 

-0.1899 

- 

-0.187 

torus 

-0.18473(4) 

- 

-0.1840(1) 

torus 

-0.183 ±0.003 

- 

-0.182(3) 

- 

-0.182 

torus 

-0.18029 


IV. SPIN-SPIN CORRELATIONS 

In this section, we examine the spin-spin correlation 
functions. Both the short-range and long-range behav¬ 
ior gives detailed information on the phases and phase 
boundaries.Since there is no long-range magnetism (ex¬ 
cept in the ferromagnetic phase, where the order param¬ 
eter commutes with the Hamiltonian), the correlation 
function is simply 

0.(z,j) = (S.-S,) , (6) 

where i and j are the indices specifying spin vertices in 
the lattice, Fig. 

A. Short-range Correlations 

To identify the bulk properties of the ground-state, we 
plot six reference correlation functions in Fig. There 
are the short-range correlations calculated for the cen¬ 
tral few sites of the 30 x 3 cylinder. The edges of 
the lattice show non-negligible boundary effects, however 
away from the boundary, the bulk correlations appear to 
be well representative of the thermodynamic limit and 
agree closely with correlators calculated using iDMRG. 
In Fig. brown stripes indicate the phase transitions 
that we have identified. 


B. Long-range Correlations 

We now consider the long-range behavior of the spin- 
spin correlators Eq. §)■ One can choose different paths 
to study distant correlators according to the lattice ge¬ 
ometry, but at long distances the spin-spin correlators 



0 

e 


FIG. 7. (Color online) Short-ranged spin-spin correlation 
functions for the ground-state on a 30x3 cylinder. Sn repre¬ 
sents the spin operator for the site n according to the lattice 
numbering in Fig. Dashed red lines are the limits of the 
correlation functions for spin-1/2 particles, —3/4 < Os{i,j) < 
1/4. Brown stripes indicate phase transitions. 

are independent of the choice of path. The Fig. shows 
correlators calculated for the path ACA as shown in the 
inset. We also calculated the correlation functions for 
a number of different paths. Up to trivial differences 
caused by the order in which different sublattices are 
listed, the results are insensitive to the path followed. 
The results suggest that the 120° and columnar states 
are quasi-long-range and the Majumdar-Ghosh state con¬ 
tains only short-ranged spin-spin correlations. 

Using iDMRG, one can directly extract the correla¬ 
tion length from spectrum of the transfer matrix. If A is 
the largest magnitude eigenvalue in the transfer matrix 






























7 




FIG. 8. (Color online) Spin-spin correlators for the THM on 
a 30x3 cylinder for trial path ACA. R{i, io) is spatial distance 
between vertices i and io in units of the lattice spacing. Two 
possible paths are shown in the inset. 

smaller than 1, then the correlation length rj is obtained 
from 

|A| = , (7) 

where oq is the size of the iDMRG unit-cell. Upon in¬ 
creasing the number of states, m, the observation of 
power-law growth of the correlation length indicates a 
gapless phase , wh ereas saturation of ?7 is a sign of a 
gapped phase.E^l^ The result for the correlation length 
of the 120°, columnar, and Majumdar-Ghosh phases are 
shown in Fig.J^ This is consistent with finite DMRG 
results of Fig. where the 120° and columnar states are 
quasi-long-range and the Majumdar-Ghosh phase has a 
finite correlation length. 


V. CHIRALITY 


In pa st decades , there has been much 
speculatiorfi^IMIZIIIIlll about the existence of chiral 
order in the 2D model. A proper chiral order param¬ 
eter will detect breaking of P and T-symmetry of the 
wavefunction while the system preserves PT-symmetry. 
This can be done by looking at order parameters or 
correlation functions that are not symmetric under P or 
T. 


We studied the chirality using two chiral order param¬ 
eters introduced below, Eq. 0 and Eq. ( [I4| ), which we 
evaluated using finite DMRG. The results are presented 
in Eig.[T0|and Eig. 0 These results show that there is no 
long-range chiral order. We also directly measured the 
parity and time-reversal symmetry of infinite size 3-leg 
cylinders, using infinite DMRG. The procedure for this 
is to calculate the overlap per unit cell of the iDMRG 
wavefunction with its conjugate or parity-reflected ver¬ 
sion. Since iDMRG works directly in the thermodynamic 


FIG. 9. (Color online) iDMRG results for the correlation 
length. T] is the correlation length, calculated from the next- 
leading eigenvalue of the transfer matrix, Eq. 0- The power- 
law growth ri oc m'‘ indicates gapless quantum critical excita¬ 
tions with power-law correlations, while the saturation of the 
correlation length in the Majumdar-Ghosh region indicates 
that this phase is short-range and gapped. 

limit, spontaneous breaking of discrete symmetries can 
occur, and this is a re liable way to detect P or T sym¬ 
metry breaking.l^ETlini xhe calculated overlap, /, between 
the P-transformed, and T-transformed, wavefunctions is 
of the order of 1 — / « I0“® per unit cell, showing that 
there is no broken P or T-symmetry. 


A. Vector Chirality 

To measure the local chirality, we use the cross product 
between vertex pairs in a plaquette, while keeping a hxed 
cyclic order of operators, 

C,[A,B,C] = SaxSb + SbxSc + ScxSa, (8) 

where [A, B, C] stands for a triangular plaquette com¬ 
posed of vertices from sublattice A, B, and C. Note how¬ 
ever that the magnitude of the local chirality is not a 
good order parameter, since it is easy to show that for 
any spin-1/2 system we have, 

3 45 

(S, X S, + S, X Sfc + Sfc X S,)2 = --MA fc + - , (9) 

where Mf^ ^ is the square of the local magnetization, 

= (S, + S, -h Sfe)2 . (10) 

Hence the magnitude of the cross product is directly re¬ 
lated to the local magnetization and has no connection 
to the chirality. 

The correlation function of the vector chirality, Oc, 
detects long-range chiral order, 

Oc{i,j-,k\i',j',k') = {Cc[i,j,k]-Cc[i',j',k']) . (11) 
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FIG. 10. (Color online) Vector chirality correlator, Oc, 
Eq. (111, on a 30x3 cylinder. Rp is the distance between the 
centers of the plaquettes in units of the lattice spacing. The 
“antiferrochiral” pattern in the 120° state can be explained 
by the tripartite symmetry of the lattice. Antiferrochirality 
is clearly broken in the Majumdar-Ghosh state, (inset) The 
path for which the vector chirality correlators were calculated. 
The P+ /P_ labels indicate the antiferrochiral ordering. 


DMRG results for this correlation function are shown in 
Fig. [T^ To calculate these correlators between desired 
plaqeuttes, we chose a path that has the maximum num¬ 
ber of crossings of plaquette vertices. This path is shown 
in the inset of Fig. 10 The origin plaquette is indicated 
in red. The next two plaquettes respectively have two 
and one common vertices with the origin while longer 
range plaquettes have none. The results of Fig. [T0|sug- 
gest that the 120° and NNN Majumdar-Ghosh states are 
only short-range chiral. There is a long-range “antifer¬ 
rochiral” pattern in the 120° state specified with P_|_/P_ 
notation in the inset of Fig. 10 which is consistent with 
the tripartite structure of the lattice. We calculated the 
vector chirality for all possible plaquettes, and all show 
the antiferrochirality of the 120° state. 


B. Scalar Chirality 

A commonly considered chiral order parameter for the 
THM is the triple product on a triangular plaquette, 
known as the scalar chirality, 

= Sa.(Sb X Sc) . (12) 

The triple product operator breaks both P and T sym¬ 
metries, and would acquire different signs for different 
plaqeuttes according to their chirality. A non-zero value 
of Ct\A, B, C\ also implies that the spins are non-planar 
on that plaquette. As a result a non-chiral and pla¬ 
nar system should acquire values close to zero for this 
triple product. Some studie^^^*^ predict that the THM 
should be chiral in some circumstances (e.g. considering 


FIG. 11 . (Golor online) Scalar chirality correlator, Ot, 
Eq. (141, on a 30x3 cylinder. Rp is the distance between 
center of plaquettes in units of the lattice spacing. The path 
chosen to calculate these correlators is the same as the inset 
of Eig. The rapid reduction of Ot to zero at long range 
indicates that the phases are non-chiral and planar. 


couplings high er than tw o-body exchange interactions), 
while the otherstrongly suggest that the quan¬ 
tum fluctuations always select a planar spin arrangement, 
so there is no chiral symmetry breaking. 

It is important to note that, the squared of the triple 
product, Eq. on a single triangular plaquette is not 
a good order parameter to measure chirality, because it 
can be shown,!^ similarly to the cross product, that for 
any spin-1/2 system, 


(S,.(S,xS,))^ = -^M2^._, + g. (13) 

As a result {Cf [A, B, C]) on a plaquette is directly related 
to the local magnetization and so is always non-zero, and 
gives no indication of the chirality. 

A diagnostic for the chirality is the correlator of the 
scalar chirality, 

= {Ct[i,j,k]Ct[i',f,k']) . (14) 

The results for the scalar chirality correlator are pre¬ 
sented in Fig. mi The path here is same as the inset 
of Fig. All phases other than FM show some short- 
range chiral correlations. However the rapid drop of Ot 
to zero at long distance is a clear sign that all phases are 
non-chiral and planar. 


VI. 120 ° ORDER PARAMETER 

Consider the classical 120 ° spin state on the triangu¬ 
lar lattice. This would be a many-body state such that 
every NNN bond is aligned ferromagnetically, while NN 
sites form AFM bonds with uniform expectation values. 
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FIG. 12. (Color online) Squared sublattice magnetization of 
120° state on a 30x3 cylinder, . This order parameter 

correctly identifies the 120° phase, consistent with the spin- 
spin correlations. Fig. 


(Si.Sj) = —1/8, throughout the lattice. This state ap¬ 
pears in a semi-classical analysis.^ 

The quantum analog of this classical 120° state can 
be constructed by positioning three spins at 120° an¬ 
gles on the Bloch sphere, forming a product state with 
long range order at wave vector Q = (47r/-\/3,47r/3), 
where the factor 2 arises from the rotation properties 
of spin-1/2 systems. The spin correlations are (Si • 
Sj) = SiSj cos 120° = —1/8 for each pair, which coin¬ 
cides with the classical value, as does the triple prod¬ 
uct Sy 4 • (Ss X Sc) = 0. The plaquette magnetization 
(Sa -I- Sb -I- Sc)^ = 3/2 is inherently non-classical. 

A suitable order parameter to detect this state is the 
squared sublattice magnetization of 120° state and can 
be constructed as, 


E 


120 ° 





47r 


47r 

X cos 

-^{x,-xe) 

COS 

- Vt') 


( 15 ) 


where No = N{N + 4)/8 is a normalization factor. 
will detect any state close to conventional 120° order. 
DMRG results for the squared sublattice magnetization 
of 120° state on a 30x3 cylinder, are 

shown in Fig. 1121 The region with non-zero values of 


^2 is consistent with the 120° phase region 


of Fig. showing that this is a good order parameter for 
the 120° phase. The value of sublattice magnetization for 
the NN model, = 0)) ^ 49% of the classical 

value, is comparable to previous calculations on the 2D 
model, 50% by ED,ESl and 40% by CCM.I^ni 

The 120° order parameter is close to maximal in the 
limit 9 —>■ —7r/2, where the ground-state tends toward the 
quantum counterpart of the classical 120° state. This 


FIG. 13. (Color online) Difference between energy per bond 
of partially polarized states and non-polarized state on a 30x3 
cylinder at 6^ = 153 °, near the critical point 8c = 152.0 °. S is 
the total spin, which is a good quantum number. The curve 
has minimum at a non-zero polarization, which indicates that 
there is a second-order phase transition close to this point. 

limit can be understood as fully saturated ferromag¬ 
netism on each sublattice via the large negative J 2 , and 
a small positive Ji induces the ideal 120° state. 


VII. PHASE TRANSITIONS AND CRITICAL 
POINTS 

In this section we pinpoint the location of the phase 
transitions and their nature, determined from the mag¬ 
netization, order parameters and spin gaps. 

The point 9 = —^ marks rapid changes in many ob¬ 
servables, consistent with a first order transition. Indeed, 
since Ji = 0 at this point, the ground-state consists of 
three uncoupled sublattices, with ferromagnetic bonds 
within each sublattice, as discussed in Sec. |III[ Therefore 
the ground-state is V/2-fold degenerate, and hence the 
120° state and the fully polarized ferromagnet coexist. 

The elementary excitations in the Majumdar-Ghosh 
chain are pairs of spin-1/2 solitons.^^ In the NNN 
Majumdar-Ghosh phase of the 3-leg ladder, the soli- 
tons in each sublattice are pinned to each other, form¬ 
ing a dislocation line. Hence the elementary excitations 
are a pairs of dislocations, with total spin S = 3. To 
study the nature of the phase transition at 9,, = 152.0° 
(NNN Majumdar-Ghosh to ferromagnet), we calculated 
the lowest-energy state in every possible total spin sec¬ 
tor. At points near the transition, we found a partially 
polarized ground-state. For example, at 9 = 153°, shown 
in Fig. the ground-state for a 30x 3 cylinder has to¬ 
tal spin S' = 3. This indicates a second-order transi¬ 
tion. We also calculated the ground-state magnetization 
around the critical point, which is shown in Fig. |14[ The 
obtained transition point indicated by the brown stripe 
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FIG. 14. (Color online) Magnetization (ground-state spin po¬ 
larization) around the critical point 6c = 152.0° on a 30x3 
cylinder. The onset of the phase transition is indicated with 
the brown stripe. 


is consistent with the correlation function results from 
Fig.0 

Accurately locating of the phase transition from the 
Majumdar-Ghosh state into the gapless columnar state 
is more difficult. Deep in the columnar phase, finite- 
size scaling of the spin gap is consistent with zero gap, 
as expected. But the finite size scaling is difficult to 
perform near the phase boundary because the finite-size 
corrections in the two phases scale differently. Hence the 
spin gap has fairly large error bars in this region, and 
the exact transition is difficult to identify. One can use 
instead the dimer order parameter, defined for this model 
as, 

= ^ E (-l)'S.-S,, (16) 


where the sum is over all NNN spins in one sublat¬ 
tice. However the dimer order parameter also contains 
large finite-size corrections. A standard procedure (al¬ 
though not common in DMRG calculations) is to use 
higher moments of the order parameter to cancel out 
low-order finite-size effects, for example using the Binder 
cumulant,!^ 

T) 


The spin gap and Binder cumulant of the dimer order pa¬ 
rameter are shown in Fig. 15 To obtain the spin gap, we 
firstly calculated the gap between S = 0 and S' = 1 total 
spin sectors for finite-length cylinders. The gap was ex¬ 
trapolated to the thermodynamic limit. Fig. 15 ^ using 
the scaling identified by Neuberger et al.,1^ which 

produces a good fit except very cfose to the transition to 
the gapped Majumdar-Ghosh phase. The Binder cumu¬ 
lant, Fig. (Mb), shows the expected behavior, whereby 




• ^ 
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5 0 

c 

^ ^ it 

/ 

/ 


/ 

/ 

/ 

/ 

3- 

1 

l- 

^ Columnar 

State 

1 

1 

^Majumdar - 
1 -Ghosh 

State 

t . ;(a) 


30 40 50 60 70 80 90 

0(degrees) 



FIG. 15. (Color online) (a) Spin gap extrapolated to the 
thermodynamic limit, (b) Binder cumulant of the dimer order 
parameter, Eq. (161, for the columnar to Majumdar-Ghosh 
phase transition. 


the value of the Binder cumulant at the phase transi¬ 
tion is independent of the lattice size (up to higher order 
corrections). The curves for 12 x 3, 24 x 3 and 30 x 3 
intersect quite closely, indicating that the transition is in 
the vicinity of 0 ~ 70.0°. 

The columnar and 120° phases are both gapless, but we 
identify the location of the phase transition from the van¬ 
ishing of the short-range order parameter shown in 

Fig. [T^ giving the transition point as 0 ~ 6.5°. 


VIII. CONCLUSION 

We have performed a comprehensive study of the 
phase diagram of the triangular J 1 -J 2 model on a 3-leg 
cylinder, using both finite DMRG and iDMRG meth¬ 
ods. There are four distinct phases in this model. All 
phases are non-chiral and planar. The 120° and colum¬ 
nar phases are gapless with quasi-long-range correlations. 
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For large J 2 > 0, the geometry of the ladder results in 
a Majumdar-Ghosh-like model with long-range dimer or¬ 
der and a two-fold degenerate ground-state. This phase 
is an effect of the restricted geometry, and only exists for 
L X 3 and L x 4 cylinders. 

Because we use a finite-width chain, the absence of 
SU{2) symmetry-breaking magnetic ordering means that 
the long-range physics is rather different to the 2D model. 
In the true 2D model, both the 120° and columnar phases 
are expected to be S’!/(2)-broken long-range ordered. 
Thus on increasing the width of the cylinder, we expect 
that the correlations will increase in magnitude and the 
gapless modes arising from the ID criticality will evolve 
into Goldstone modes associated with the broken sym¬ 
metry of the order parameter. 

The short-range physics and structure of the phase di¬ 
agram of the 3-leg ladder agrees closely with known re¬ 
sults for the 2D model, especially in the small J 2 region. 
We find a transition from 120° to columnar phases at 
9c — 6.5°, close to the classical value. Further studies 


on larger width cylinders have clarified that between the 
120° and columnar state there is a spin liquid region,!^ 
consistent with t he re cent results of quantum Monte 

Carlo calculations.l^^lll 

The boundary between ID and 2D physics in this 
model is rather rich, and this suggests that the physics 
arising from restricting geometry to finite-width ladders 
presents a fruitful direction for future investigation, and 
may explain some novel properties of molecular solids.l^ 
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